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Abstract 

Hereby we complete the proof of integrability of the Lax systems, based on 
pseudo-Riemannian coset manifolds G/H*, we recently presented in a previous paper 
arXiv:0903.2559 . Supergravity spherically symmetric black hole solutions have 
been shown to correspond to geodesies in such manifolds and, in our previous paper, 
we presented the proof of Liouville integrability of such differential systems, their 
integration algorithm and we also discussed the orbit structure of their moduli 
space in terms of conserved hamiltonians. There is a singular cuspidal locus in 
this moduli space which needs a separate construction. This locus contains the 
orbits of Nilpotent Lax operators corresponding to extremal Black Holes. Here we 
intrinsically characterize such a locus in terms of the hamiltonians and we present 
the complete integration algorithm for the Nilpotent Lax operators. The algorithm 
is finite, requires no limit procedure and it is solely defined in terms of the initial 
data. For the SL(3, M)/S0(1, 2) coset we give an exhaustive classification of all 
orbits, regular and singular, so providing general solutions for this case. Finally we 
show that our integration algorithm can be generalized to generic non-diagonalizable 
(in particular nilpotent) Lax matrices not necessarily associated with symmetric 
spaces. 
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Landau Program. 



1 Introduction 



As announced in the abstract, this paper completes the results of our very recent paper 
[1] where we discussed the classification of orbits for G/H* Lax systems, using the con- 
served hamiltonians arising from the underlying Poisson structure and adapted Kodama's 
integration algorithm [2] to this case. For the motivations, the general setup and a more 
extended bibliography we refer the reader to that publication. Assuming the framework of 
[T] we just focus on the problem of nilpotent orbits, that have a distinct physical interest 
since they correspond to extremal Black Holes. The authors of [3j , in parallel to our paper 
[Ij, have also independently adapted Kodama integration algorithm to the treatment of 
G/H* for the case of non nilpotent initial conditions and have presented some particular 
solutions. 

The results presented here are three: 

A. The classification of exceptional orbits of Nilpotent Lax operators. 

B. The explicit integration algorithm for such operators. 

C. A generalized formula giving solely in terms of the initial data Lq the integration of 

Lax equation for arbitrary non-diagonalizable and diagonalizable Lax matrices L, 
not necessarily representing symmetric spaces. This is specially useful if Lq is non- 
diagonalizable (in particular nilpotent) since this is the only available integration 
formula. This result is presented in the appendix. 

Although the presented integration algorithm is general, to illustrate it, we heavily rely on 
the simple SL(3,M)/S0(1, 2) example, already employed in [1]. The same paradigmatic 
example is utilized to explain the classification of orbits. 

2 The singular orbits of Nilpotent operators 

The integration algorithm determining the solution of the first order differential equations 
for the tangent vector to a geodesic in a manifold G/H* was discussed in |lj where their 
Liouville integrability was proved through the construction of the preserved hamiltonians, 
constructed with the following procedure. Given the Lax operator L, fulfilling the 1]- 
symmetry condition 

{Lvf = Lv (2.1) 

with 

r/ = diag(-l,+l,...,-l,+l,+l,... + l) , p<q ; p + g = N (2.2) 

the complete set of pn functions i)a that are involutive with respect to the Lie-Poisson 
bracket defined on the solvable Borel Lie algebra of N x N upper triangular matrices is 
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enumerated by an ordered pair of indices 



a 



{a,b) 



(2.3) 



where: 



a 





b = l,...,N-2a. 
The functions i)ab can be iteratively derived from the following relation: 

det {{L - X)ij : a + l<z<N, l<j<N-a} 

N-2a \ 



(2.4) 
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N-2a-b 



0, 



(2.5) 



where, by definition, Sao is the coefficient of the power A^~^". Among the conserved 
hamiltonians there are polynomial ones that depend only on the eigenvalues of the Lax 
operator and rational ones that instead depend also on the initial twisting data (the H*- 
rotation of the Lax operator from its diagonal form in the case of real eigenvalues or its 
block-diagonal normal form in the case of complex eigenvalues). There are also among 
the i)ab Casimir functions that have vanishing Poisson brackets with all the canonical 
variables. 

As we extensively discussed in the classification of spectral types, which corresponds 
with the classification of normal forms given in can be mapped into the foliation of 
the space of Lax orbits spanned by the conserved hamiltonians of polynomial type. 

For instance in the case of SL(3,M)/S0(1, 2), where the general form of the Lax 
operator is given by [1] 



leading to the following system of differential equations: 
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there are just a quadratic and a cubic hamiltonian 

, (2.8) 

\yiitf + \Y,{tf - -^,{tf + \Y,{tf - \Y,{tf , (2.9) 

Y^jt)^ Y,{t)'Y2{t) ^ YsjtfY.jt) ^ Y.jtfY^jt) Y.jtfY^jt) 
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while the rational hamiltonian is the following one: 

The space of orbits is separated in two distinct regions by the value of the discriminant 

A = -12i)l + 81i)l. (2.12) 

• In the region where A < we have three distinct real eigenvalues. 

• In the region where A > there is one real eigenvalue and a pair of complex 
conjugate eigenvalues. 

• The locus A = corresponds to orbits admitting an enhanced symmetry, except at 
the cusp. 



2.1 The orbits of the J = 2 representation of S0(l,2) 

The easiest way to discuss the space of orbits for the SL(3, ]R)/S0(1, 2) case is to recall that 
the coset generators span the five dimensional J = 2 representation of S0(l,2), drawing 
an analogy to the well known case of orbits of the J = 1 representation. In the latter case 
we deal with vectors and the space of orbits decomposes into three sectors: 

1. The orbits of time-like vectors {v , v) > which admit SO (2) as stability subgroup 
and therefore span the coset S0(1, 2)/SO(2). 

2. The orbits of space-like vectors {v , v) < which admit S0(1, 1) as stability sub- 
group and therefore span the coset S0(1, 2)/S0(l, 1). 

3. The orbits of light-like vectors {v , v) =0 which admit no stability subgroup and lie 
on the light-cone which separates the two regions of time-like and space-like orbits. 
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Figure 1: Orbit structure in the [)2, i)3 plane. The cuspidal point corresponds to the 
Nilpotent orbits 

In the J = 2 case the role of the light-cone is played by the A = locus of vanishing 
discriminant. The orbits with A > have the eigenvalue structure we have just discussed 
and admit no stability subgroup. This can be directly verified using the generators of 
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S0{1, 2) given in [T] for the J = 2 representation and recalled here for convenience: 
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(2.13) 

In this way the original five dimensional space is parameterized for A > orbits by the 
values of the two hamiltonians f)2 , J) 3, which label the orbit, and by the three parameters 
of S0(l,2) that span it. Similar conclusion one draws for the A < 0. Apart from the 
difference in spectral type these orbits behave exactly in the same way. There is no 
stability subgroup. 

On the other hand the locus A = contains 5-vectors that always admit a stabilizer 
S, either belonging to the compact S e so (2) subalgebra or to a non-compact subalgebra 
S G 50(1, 1) of so(l,2). For these orbits the discriminant A vanishes, yet the individual 
hamiltonians f)2 and f)3 are generically different from zero. In this way each of the A = 
orbits spans either the S0(1, 2)/SO(2) or the S0(1, 2)/S0(l, 1) coset. 

There is finally the cuspidal locus displayed in fig jl] where both [)2 and f)3 are zero. The 
five-vectors lying in this locus correspond to Nilpotent non-diagonalizable 3x3 matrices 
as, for example, the following: 



/o 

1 



-1 o\ 

1 

1 Oy 



(2.14) 



This provides an intrinsic characterization of the Nilpotent orbit: vanishing of both poly- 
nomial hamiltonians. 

Another interesting and equivalent algebraic characterization is the following. The 
2-dimensional space of non-compact generators of the Lie algebra so(l,2) given by the 
span of Jl and J3 contains operators whose eigenvalues are necessarily: 



-2i, -i,0, +i, +2> 



(2.15) 
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where ^ is an arbitrary normalization of the operator. All the Nilpotent operators L 
belonging to the singular cuspidal orbit Cq are eigenvectors of one of these non-compact 
generators G with a finite eigenvalue /x 7^ 0: 



On the contrary one can easily verify that no five-vector corresponding to regular orbits 
can be eigenvector of any of the so(l, 2) generators. 

In view of these remarks we see that the A = locus is also alternatively characterized 
as the set of operators that are eigenstates of some so(l, 2) generator. Null eigenstates oc- 
cur along the two branches and give orbits with enhanced symmetry (stability subgroup). 
Eigenstates of non- vanishing eigenvalue occur only at the cusp. 

Conjecture 2.1 For higher groups with more hamiltonians, the secular equation has typ- 
ically degree higher than four and the discussion of discriminants becomes unavailable. 
Yet the second algebraic characterization of cuspidal orbits of Nilpotent operators that, by 
definition, have all vanishing eigenvalues and are not diagonalizable, remains viable and 
appears very promising. Indeed we conjecture that the orbits of such Lax operators can be 
found by this method, namely searching for eigenstates of the non-compact generators of 



We stress that in the case of supergravity billiards [5], [6] the isotropy group H is 
compact and therefore there are no real Lax operators that can be eigenstates of any 
generator of H. It is precisely the pseudo-Riemannian nature of the coset G/H* what 
allows for the appearance of nilpotent Laxes. 

For this singular, isolated class of initial conditions the integration algorithm was not 
provided in [1] since one of its ingredients consists of the diagonalization of the initial 
Lax. The next section fills this gap showing that the solution can be directly and simply 
constructed in terms of the initial Lax Lq. 

3 The integration algorithm for Nilpotent Lax oper- 
ators 

Let us assume that the initial data for the Lax equation 




(2.16) 



d_ 
dt 



L{t) + [L 



> 



L<,m] = 




are provided by an r^-symmetric matrix: 



L(0) = U 



(3.2) 
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which is also nilpotent and therefore not diagonahzable. Let us call n the minimal degree 
of nilpotencj|^ 

(3.3) 



Lo = 0. 



The complete solution for L{t) that admits such an initial condition is constructed in the 
following way. First define the following building blocks: 



1. The following N x N matrix function: 



71—1 _ 



(3.4) 



*;=0 



2. The minor functions constructed in the following way: 

/ Ci,i(t) ... Ci,,„i(t) \ 



miit) = 1 



Ck,i(t) ... Ck,i-i{t) 



V Ci,i(t) ... Ci,,_i(t) / 



, 1 < k <i; 2 <i < N , 



(3.5) 



where the hats on the entries corresponding to the fc-th row mean that such a row 
has been suppressed giving rise to a squared {i — 1) x [i — 1) matrix of which one 
can calculate the determinant. 



3. The determinant functions: 



Di(t) = Det 



/ Ci,i(t) ... Ci,{t) \ 
\ C,,i(t) ... a,{t) J 



Do(t) :=1. 



(3.6) 



Then in terms of them we write the entries of the izme-evolving Lax operator in the 
following way: 



Lpq{t) 



P 1 



S)p(t)Dp_i(t)D,(t)D,.i(t) 



^In the SL(3,]R) case the degree of nilpotency is either n — 2 or n = 3. It seems that the same occurs 
also in the case of more general cosets G/H* occurring in supergravity theories. However the algorithm 
we present applies to any finite degree of nilpotency n. 
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As an example of the method we give the exphcit form of the Lax operator which is 
produced by the integration with initial conditions provided by the operator Q displayed 



in (2.14): 



( 



\ 



2t 



l-2f2 
l-2t2 







Vl+2t2 
1-2*2 
4t_ 
1 



4*4- 

Vi- 



1+2*2 





Vl-2t2 
1+2*2 

2t 
1+2*2 



(3.8) 



This is a representative of the orbit where the Lax operator has degree of nilpotency 
n = 3. It corresponds to the eigenvalue /i = 1 of the non compact 5o(l,2) generator 

R(J3). 

It is quite interesting to employ the algorithm for the square of the matrix Vt which 
is an eigenvector corresponding to eigenvalue /x = 2 and has degree of nilpotency n = 2. 
The result is a very simple solution with a completely different analytic behaviour 



L2{t) 



l + 2t 



-1 



1 



-l\ 



1 / 



(3.9) 



3.1 Summarizing 

We conclude that the two above presented solutions are the representatives of the two 
nilpotent inequivalent orbits corresponding to degree of nilpotency n = 3 and n = 2, 
respectively. Any other solution of the same classes is obtained by applying our algorithm 
to an S0(1, 2) rotation of the corresponding initial data. All in all the three-dimensional 



manifold of nilpotent Lax operators is provided by those matrices (2.6) whose five pa- 
rameters Yi satisfy the two constraints 1)2 (^) = f)3(^) = 0, according to the definitions 



(2.9-2.10) of the quadratic and cubic hamiltonians. 



The other exceptional orbits are arranged along the two branches of the A = locus 
and each of them spans either an S0(1, 2)/SO(2) or an S0(1, 2)/S0(l, 1) coset manifold 
since there is either an so(l, 1) or an so(2) stability subalgebra. These orbits correspond to 
diagonalizable Lax operators and are covered by the algorithm presented in our previous 
paper [1]. 



4 Conclusions 

In this note we have shown how the integration algorithm for Lax equations can be ex- 
tended also to the case of non-diagonalizable (see, also the appendix) initial data. Appli- 
cations to the construction of extremal BPS solutions is postponed to future publications. 

Aknowledgments We would like to express our gratitude to our frequent collaborator 
and excellent friend Mario Trigiante for the exchange of useful information we had with 
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him in these last few days. He pointed out to us the relevance of the nilpotent case and 
formulated the problem that we were able to solve in complete generality. 

Note added in revised version We recall that a particular solution corresponding to 
a specific nilpotent initial condition for the SL(3, ]R)/S0(1, 2) case was presented in the 
revised version of ^ which appeared on the hep-th ArXiv the very same day as the first 
version of the present paper, where the problem obtained its general solution. 



A A generalized formula 

Although the physical or geometrical applications of such an equation are not immediately 
evident to us we remark that one could consider a generalization of the Lax equation for 



symmetric spaces G/H by writing the very same differential condition (3.1) imposed on 



a matrix L{t) which is not required to satisfy the ?7-symmetry conditions (3.2). In other 



words we can extend Lax equation to matrices L(t) that are not necessarily elements of 
the orthogonal complement K to a subalgebra H* C G C s[(N,]R). 

It is interesting that we can integrate Lax equation in a completely general way with 
arbitrary initial Lax matrices Lq. This is done by substituting eq.(3.7) with the following 



one: 



k=i e=i 



where 



m,k{t) := (-l)*+'=Det 



Ci,fc(t) 



l<k<i ; 2<i<N ; ^ii(t) = 1 , 



(A.2) 



and all other definitions remaining the same as in the main text. In eq.(A.2) the hatted 



k-th. column is deleted just as in eq.(3.5) it was deleted the k-th row. 



This expression gives the solution to Lax equation (3.1) for the case of generic initial 
conditions. It applies in particular to the case of diagonalizable initial matrices Lq. In 



that case, however, eq.(A.l ) is only formal and not too useful since the function C{t) (3.4) 
is constructed by a matrix exponentiation and involves summing an infinite series. On 



the contrary for nilpotent matrices Lq eq.(A.l) provides a useful explicit general integral. 



The class of nilpotent matrices is much larger than the class of nilpotent //-symmetric 



matrices (3.2), the latter being associated with symmetric spaces. 



Although the physical meaning of Lax equation in the more general context of nilpo- 
tent, but not necessarily //-symmetric matrices (3.2) is so far unknown, yet it is worth 
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mentioning the existence of this generahzed integration formula whose apphcations will 
certainly be discovered soon. 

Now, let us derive another equivalent representation of equations (A.l) which will 
be useful in case of both nilpotent and non-nilpotent, ?7-symmetric and ?7-non-symmetric 
initial data Lg entering these equations. 

At first, we recall that both diagonalizable and non-diagonalizable (in particular nilpo- 
tent) generic initial conditions for the Lax matrix Lq can be represented in the Jordan 
normal form 

Lo = ^{0)J^{Oy' (A.3) 

where $(0) is an invertible N x N matrix, JT" is a block-diagonal N x N matrix with the 
da X da matrix sub-blocks 



J 



f Jx, ... \ 

■•• ■•• : 

; ■•. ■•. 

\ ... Ja„ / 



1 







... 

■•• 

■•. 1 

A„ 



\ 



(A.4) 



a 



l,...,m, l<m<A^, (ia>l, d^ 



N 



a=l 



and Xa are constants. Substituting eq. (A.3) into equations (A.l) and (3.4), the latter 
can identically be represented in the following equivalent form: 



L{t) = 

C{t) = $(0)6-2*^ $(0)-^ 



(A.5) 
(A.6) 



where 



N 



I C,^r{t) ... Ci,_i(t) $1,(0) \ 



5^Det 



S3 



\ Ci,i{t) ... C,,,_i(t) $,,(0) J 



( 



N 



i^-'),. it) 



v/2'^(t)S),_i(t) ^ 



\ 



Ci-l,l(t) . . . Ci_i^j(t) 

V (<^-^).i(o) ... ($-^.(0) ; 



(A.7) 



and the functions Di{t) are defined in eq. (3.6). 
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For the particular case of diagonalizable initial conditions, i.e. when m = N {da = 1, 
N), the corresponding matrix J' (A. 4) becomes a diagonal matrix of the A^- 
1, . . . , A^), and equations (A. 5- A. 7) reproduce the general solutions 



[a 



a = 1, . 
eigenvalues A 

to Lax equations (3.1) for the case of generic diagonalizable Lax matrices derived in [7| 
for the first time. What concerns expressions (A.5-A.7) in the case m N, they are 



a generalization of the above-mentioned general solutions to the case of generic non- 
diagonalizahle Lax matrices constructed, to our best knowledge, for the first time in the 
present paper. 



When deriving eqs. (A.5-A.7) starting from eq. (A.l), we did not use any properties 



of the Lax operators Lq but its rather general representation (A.3). In case of diagonal 



matrices J we proved the correctness of derived expressions (A.5-A.7) by their estab- 



lished relationship with the earlier known solution [7] of Lax equations (3.1), as it was 



cosets at 



already mentioned above. We also verified their correctness in case of ^^^^ 
A^ = 2 , 3 , 4 and 5 for non-diagonalizable nilpotent matrices J . Altogether, this gives an 
evidence in favour of correctness of eqs. (A.5-A.7) in case of generic matrices J (A.3) as 
well. 



The constructed general solutions (A. 5 -A. 7) are explicitly parameterized by the ini- 
tial data encoded in the matrices J and $(0). In order these solutions could be ex- 
pressed in terms of elementary functions, the exponential e~*"^, entering into relations 
(A. 6) and(A.7), has to admit a closed form in terms of elementary functions. Obviously, 
it is a very simple task for diagonal and nilpotent matrices J . Thus, in the latter case 



all eigenvalues are zero since J in (A.4) is a nilpotent matrix, and the expansion of the 
exponential e~*"^ in power series with respect to J terminates at some finite order. It 
is more complicated but in principle a solvable task for many interesting, more generic 
non-nilpotent and non-diagonalizable cases. 



Thus, for an example in case of ?7-symmetric Lax operators Lq (3.2), is a pseudo- 
orthogonal matrix, G SO{p,q), i.e. = ri^'^{t)ri. In this case the non- 
diagonalizable matrices J7 have a very simple structure [1] : 



J = Q + Nil, [Q, Nil] = 



where Nil is a nilpotent matrix and Q is a block-diagonal matrix with 2x2 and 1x1 
sub-blocks [HII]. Therefore, the exponential e~*"^ can be factorized 



..Nil 



(A.9) 



and both exponentials, entering on the r.h.s. of the latter relation, can easily be expressed 
in terms of elementary functions. A detailed discussion of this interesting case will be 
presented elsewhere. 
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